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Abstract. We extend an above barrier analysis made with the Schrodinger equation to 
the Dirac equation. We demonstrate the perfect agreement between the barrier results 
and back to back steps. This implies the existence of multiple (indeed infinite) reflected 
and transmitted wave packets. These packets may be well separated in space or partially 
overlap. In the latter case interference effects can occur. For the extreme case of total 
overlap we encounter resonances. The conditions under which resonance phenomena can 
be observed is discussed and illustrated by numerical calculations. 



PACS. 03.65.Pm - 03.65.Xp 



I. INTRODUCTION 

This paper considers the above barrier solutions of the Dirac equation for an electrostatic one dimen- 
sional (z axis) potential, 

!0 , z < (region I) , 
V , < z < I (region II) , 
, z > I (region III) . 

It is very difficult and probably even confusing to treat, in a single article, all interactions of plane 
waves or wave packets with a barrier potential using the Dirac equation. This is because the physical 
content depends upon the energy of the incoming (particle) wave. In the figure below we depict the 
three potential regions. Also shown is the energy divided into three zones. 
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The upper energy zone, E > V — to, is that of interest to this work and involves diffusion phenomena. 
In the so called Klein zone [1], E < V — to, oscillatory solutions exist in the barrier region. These 
are antiparticles [2-6] . Indeed, antiparticles see an opposite electrostatic potential to that seen by the 
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particles and hence they will see a well potential where the particles see a barrier. The antiparticles 
thus live above the well potential and are legitimately oscillatory in form. In the tunnelling zone only 
evanescent waves exist in the barrier region [7,8]. Of particular interest here is the possibility of an 
Hartman-likc effect [9-11]. 

No barrier analysis can be interpreted without an understanding of the step potential results. In 
a latter Section, we shall argue, for the step, that while pair production is innate in the Klein zone 
there is no pair production for above barrier diffusion where the non-relativistic limit reproduces the 
standard Schrodingcr solutions. There is of course no Schrodingcr limit for the Klein zone. The Klein 
zone of the step is characterized by the Klein paradox in which a reflection probability greater than 
the incoming probability occurs. The excess particle number (or charge) implies pair production. For 
above step diffusion this paradox does not exist and hence pair production is absent. 

For the case of above barrier diffusion, we shall demonstrate in detail the equivalence of a two- 
step calculation and the barrier results. Such an equivalence has previously been demonstrated for 
the Schrodinger equation [12]. This method of calculation which employs the simple step results is 
not new. Multiple step calculations have, for example, been used in the WKB applications to various 
potentials [13] and for the barrier potential in the Dirac equation in the Klein zone [14]. However, these 
authors employ the method only as a convenient mathematical tool. We, on the contrary, emphasize 
its physical significance in terms of multiple reflected and transmitted peaks (when the barrier width 
is much greater than the wave packets widths). The exit times for each can be calculated using the 
stationary phase method (SPM) [15,16]. They occur with time intervals equal to twice the barrier 
width divided by the group velocity over the barrier. 

The combined summed expressions for the reflection and transmission coefficients contain the well 
known resonance phenomena. We shall discuss under which conditions this resonance effect occurs 
and when the effect "breaks up" and the multiple peaks appear. Some numerical calculations will help 
us to illustrate this transition from an effective single outgoing wave packet (coherence) to essentially 
independent multiple wave packets (de-coherence). 

In the next Section we will give the necessary formulas and conditions assumed for our calculations. 
In Section III, we consider the plane wave solutions for a step potential and more specifically for three 
related but distinct steps. One upward step at z = and two downward steps at z = and z = I. The 
extra phases that appear in the last case are essential for the calculation of the times of the outgoing 
peaks. We will then calculate the back to back step potential obtaining the individual amplitudes 
for the (infinite) reflected and transmitted terms. Each can be associated with a wave packet after 
integrating with a suitable convolution function. In Section IV, we calculate directly the plane wave 
solutions for the reflection and transmission coefficients for a barrier and evidence their equivalence 
to the sum of the results from the previous Section. Section V discusses the question of resonance and 
points out that this phenomena requires specific conditions to occur. We conclude in Section VI with 
a resume of our results. 



II. DIRAC SOLUTIONS IN A CONSTANT POTENTIAL 

The free Dirac equation reads 

{i^df, - m)^{r,t) = (1) 

with the gamma matrices satisfying {7^, 7"} — 2g /iI/ . It has four well known independent plane 
wave solutions 

u {l - 2) (p) cxp[ip ■ r - iEt] , (E = \E\) and u {3 ' 4 \p) cxp[ip ■ r - iEt] , (E = —\E\) , 
where \E\ = ^Jp 2 + m 2 . Using the Pauli-Dirac set of gamma matrices 

-1 ) and 7= {-a 
the spinorial solutions are (for polarized states in the z-direction) 

/ X (s) \ / a ' P v (s) \ 

u^(p)=N I <r p (s) and u^ip) = N \E\+m A , (2) 

\E + m X J V X (s) J 
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where s = 1, 2 and 

X W =(J) and X (2) =(?) , 

with N = \J (\E\ + m) the covariant normalization choice such that uM t u ( s ) = u( s+2 ) t u ( s + 2 ) = 2 
We recall that the free Dirac Hamiltonian is 

= + (A -£) ■ ( 3 ) 

The different signatures of the energies between s and s + 2 spinors follows readily from the alternative 
non-covariant form of the Dirac equation 

id t V(r,t) =H *{r,t) (4) 

the left hand side yields E$>(r,t) in all four cases while the right hand size yields ± \E\ as the case 
may be. Actually, for ^^(r, t) we appear to get a simple identity because we have conventionally 
used E rather than \E\ in its spinor representation. However, when p = (the rest frame case) the 
Hamiltonian reduces to mj" and the above equation yields E — + m for s = 1, 2 and E = —m for 
s + 2 = 3,4. These solutions are oscillating solutions, valid when E > m or E < — m. There are 
also evanescent solutions obtained from the above with the substitution p — > i p for which the spatial 
dependence becomes exp[=pp] and occurs when — m < E < m. 

Since we shall need the solutions for step and barrier potentials, we rewrite the above solutions in 
the presence of a constant potential V a . Consider an electrostatic potential — (A , 0) included (via 
minimal coupling) in the Dirac equation 



or 



The spinorial solutions are now 



(j7"9,-e7%-m)*(r,i)=0 
»Wr,t) = (H„ + V a )V(r,t) 



H a V(r,t) = (E-V )#(r,t) . 



(5) 
(6) 

(7) 



where V — — e A (charge — e). For a stationary solution ^(r, t) oc exp — i Et], we obtain 



u^(q;V ) = N "I W , (8) 
\E-V + m A / 

for E — V > m, and 

/ a —3- vW\ 

u^(q;V )=N f \E-V \+m X j, (9) 

for E — V a < —m, where (E — V ) — q 2 — m 2 and N = y/ (\E — V \+ m). This shows that the two 
set of solutions and w^ 2 ) are not determined by the sign of E but by whether E > V + m or 
E < V — in. the latter being the defining condition for the Klein zone. Hence, E may be fixed but 
the solutions depend upon whether in any given region the energy is above or below the potential. 

In any constant potential region only two solutions exist for a given E (be they oscillatory or 
evanescent). This fact is essential for a standard plane wave analysis. Because the Dirac equation is 
a first order equation in the spatial derivatives, for a step-wise continuous potential only continuity 
of the wave function exists. There is no lose of information w.r.t. the Schrodingcr continuity (wave 
function and its spatial derivative) since continuity of the " small" component of the Dirac spinor yields 
in the non-relativistic limit the continuity of the Schrodinger wave function's spatial derivative. Thus, 
continuity of the Dirac wave function implies four conditions at each interface, one for each spinor 
component. However, since the sign of the momenta are a priori arbitrary, we have two conditions 
which determine the spinor. Hence, for consistency, there can be only two independent spinors in each 
region and these must correspond to different helicity states or orthogonal combinations of them. 
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III. THE TWO STEP APPROACH 



Let us treat the barrier diffusions as an application of a two step process. One upward step at z — 
and one downward step at z = I. We assume an incoming above-potential, positive helicity, plane 
wave state given by (the normalization of the spinor is inessential here) 



1 



P 



\ 



E + m 




exp[i (pz — Et) ] 



with E = y/p 2 + m 2 . We shall need three step diffusions, the first at z = both for incoming waves 
from the left and from the right. The second one at z = I only for incoming waves from the left. For 
reflection and transmission from an upward/downward potential we shall use the suffixes +/— ■ 
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For the first step, the solutions in region I and II are respectively given by 
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q 



exp[ip z 
\ 



R + (0) 



1 



V 



\ 



E + m 




exp[ — ipz] > exp[ — iEt] 



E - V + m 




(10) 



exp[igz] cxp[ — iEt] , 



with q = y/ (E — V ) 2 — m 2 . All helicity (spin) flip terms (the other independent spinor solutions) here 
and elsewhere turn out to be absent, so we exclude them a priori for simplicity. Continuity at z = 
yields 

l + i? + (0)=T + (0) and 1 - R + (0) = T + (0) q {E ± " } r = T + (0) - , 

p{E—V + m) a 

where a = \f(E — V a + m) (E — m) and b = y/ {E — V — m)(E + m) (a > b). Thus, we obtain 



Step 1 



R+(0) - (a-b)/(a + b) 
T + (0) = 2a/(a + b) , 



(11) 



from which it follows that 



l-\R + (0)\ =~\T + (0)\ . 



While \R + (0) | is obviously the reflection probability, the "transmission" probability must therefore be 

2 

(b/a) |T+(0)| . A demonstration of this for wave packets (albeit with the aid of some approximations) 
is given in the Appendix. Since |i? + (0)| < 1, we do not have a Klein paradox here so no pair production 
is involved. 

To understand physically the significance of single reflection and/or transmission coefficients one 
needs to use the stationary phase method for wave packets (normalized convolution integral of plane 
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waves). This method, introduced by Stokes and Kelvin [1], estimates the position of the maximum 
of the wave packet by using the simple concept that, far from the vanishing derivative of the phase, 
the argument of the convolution integral oscillates many times and produce destructive interference. 
Consequently, the maximum of the wave packet occurs where the derivative of the phase vanishes [16]. 
For example, consider an amplitude 

= J dkg(k) exp[i0(k)] cxp[i ( k z - E t ) ] 

modulated by a real function g(k), with a single steep maximum at k . The time-space relationship 
for the maximum (maxima if <f> has spinor components) is given by 



t = 



dk 
dE 



d0_ 

dE 



the zero suffix means the parenthesis is calculated at k = k . The delay factor is (dd/dE) . The 
group velocity is v g = (dE/dk) a . In our analysis the primary modulation function, g(p), is that of 
the incoming wave packet 



dp g(p) 



1 


P 



E + m 




cxp [ i ( p z — Et)] . 



(12) 



The "effective" modulation function in any given region is then given by g Aui where A is the plane 
wave amplitude (e.g. R+(0) or T+(0)) and it, stands for the spinor clement considered. For each 
separate spinor component one must calculate the group velocity and eventual delay times. However, 
with our choice, all non zero spinor components are real, so that there is no contribution to 8 from 
them. Thus, with real i? + (0) and T + (0) we have no time delay. It is true that the group velocities 
depends upon the spinor components momentum dependence (which shifts the value of fc ). However, 
this should be negligible for a very sharply peaked modulation function. It is demonstrably negligible 
in the two limits: Non-relativistic (NR) where one can simply ignore the small components; and in 
the ultra-relativistic (UR) when the ("small") component q/(E — V„ + m) tends to one. 

Returning to our calculations, we give without detail the results for the other two steps. For step 
2, we find 

' R_(l) = [(b-a)/(a + b)] cxp[2iql] , 
T_(0 = [2b/(a + b)] exp[i(q-p)l] , 



Step 2 



(13) 



and for step 3, 



Step 3 



R-(0) = (b-a)/(a + b) 
T_(0) = 2b/(a + b) . 



(14) 



The first of the barrier transmitted amplitudes is thus obtained by multiplying the "transmitted" 
amplitude of step 1, T + (0), by the transmission amplitude of step 2, T_(l), 

[4ab/ (a + b) 2 } exp[i(q-p)l] . 

The exit time is calculated after including the plane wave phase in region III, 



t = 

where 9 = (q — p)l. Thus, at z = I we find 

t = 



dp 
dE 



d6 
dE 



dp 
dE 



dq_ 
dE 



dp 
dE 



1 = 1^- 
1 dE 



This is just the time for a wave packet in region II with group velocity {dE /dq) a to travel a distance I 
(barrier width). At step 2 there is also a reflected amplitude given by T + (0)R_(l). The corresponding 
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wave packet travels back towards z — and the second reflected wave exits into the left region I, with 
amplitude (unmodulated) T + (0)R_(l)T_(0) at the expected time 

t = 2 {dq/dE) a I . 

The whole procedure may then be repeated ad infinitum. 

Below, we list the first few individual reflected and transmitted waves together with the expressions 
for the general n th wave. 



Ri = R+(0) 



a — b 
a + b ' 



R 2 = T + (0)R_{l)T_(0) = 4ab< t a) exp[2»gl] , 

(a + by 

R 3 = T + (0)R_(l)R_(0)R_(l)T_(0) = 4 ° 6(5 ~° ) exp[2iql] ( h —^r exp[iql] 



(a + b) 3 



a + b 



Rn = T + (0) R_ (0 [ R_ (0) R_ (I) ]"- 2 T_ (0) = 4 - 6 (& a) cxp[2iql] (^—^ cxp[iql] 



(a + b) 3 



a + b 



For the transmitted amplitude, we have 

Aab 



T 1= T + (0)T_{1) 



(a + b) 

T 2 = T + (0)R_(l)R_(0)T_(l) = 
T 3 = T + (0)[R_(l)R_(0)fT_(l) 



cxp[i(q-p)l] , 



exp[i(q-p)l] ( — 4 exp[iql] ) , 



(a + b) 2 " ' \a + b 

= 7~7T^ C MH<1-P) 1 } \ cxp[iql] 
(a + b) 2 \a + b 



(15) 



T n = T + (0)[R_(l)R_(0)] n - 1 T4l) 



Aab 

(a + b}- 



exp[i(q-p)l] 



b — a 
a + b 



exp [iqt 



(16) 



At each step one can check (in accordance with our previous discussion) that probability is conserved. 
If the individual wave packets are well separated (see the following Section) the probability of, say, the 
n th transmitted wave will be just \T n \ since it travels in a potential free region. A straightforward 
calculation then shows that the total transmission probability is, 



16 a 2 b 2 

(a + by 









1 - 




"VI 




■v J 



2ab 

a 2 + b 2 



(17) 



It is to be noted that this sum is independent of the barrier width I. A similar calculation can be 
performed for the region I reflected probabilities. This yields 



n=l 

Consequently, as expected 



a-b\ 2 16 a 2 b 2 (b -a) 



{a + by 

OO 

J2{\Rnf + \T n \ 2 }=l 





/»- 




1 - 




°V1 






v J 



(a - b) 2 



(18) 



(19) 



Our multiple peak interpretation is thus consistent with overall probability conservation. Finally, we 
observe that the time interval At between two successive outgoing peaks, in either of the potential 
free regions is 
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IV. THE BARRIER ANALYSIS 

Let us now perform the standard stationary plane wave analysis for the barrier - again neglecting a 
priori (for simplicity) spin flip. 



Region I: 



z < 0. 



P 



exp[ip z] + R 



Region II: < z < I , 



Region III: I < z , 



Continuity at z = yields 



i E + m . 

V o J 

< I ^ 

q 



p 



E + m . 


/ 



A 



T 



E - V + m. 




V 

( I \ 



exp[iq z] + B 



cxp[ — ipz\ 



1 \ 





E - V + m 




exp[ — iqz\ 



cxp[ip z] 



E + m 

V o / 



1 f a + b a-b\ ( A 



2a \ a-b a+b I V B 



From continuity at z = I, 



Consequently, 



1 1 \ / A cxp[igZ] 
1 -1 M B exp[-iql) 



1 / exp[ — iql] 



jh) T CX P[*P ; ] 



V 1 - 1 

cxp[ip I] 



lb ) T CX P[*P 



2 b 



(b + a) cxp[ — iql] 
(b — a) exp[i ql } 



T 



Using this equation to eliminate A and B from the continuity equation at z = gives 
f \ exp[ip I] 



R 



Aab 



(a + b) 2 exp[ — i q I ] — (a — b) 2 exp[ i q I ] 
(a 2 — b 2 ) cxp[ —iql] — (a 2 — b 2 ) exp[i ql] 



T 



Whence, 



and 



R = i{b 2 -a 2 ) sm[ql] / [ 2 a b cos[ql] - i (a 2 + b 2 ) sin[ql] ] 



(20) 



(21) 



T = 2ab cxp[-ipl] / [2ab cos[q/] - i(a 2 + b 2 ) sin[g/]] . 
These amplitudes satisfy, 

\R\ 2 + \T\ 2 = 1. (22) 

If treated as single wave packets (one for R and one for T) the peaks emerge at a common time 
determined by their common denominators 



in — t'i 



AO 
dE 



n a 2 + b 2 

with tanO = : — tan ql 

2ab lH J 
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Note that the phase factor exp[ — ip I] in T has been cancelled by the plane wave factor exp[ip z] 
calculated at z = I. The i factor in R is momentum independent so it does not contribute to the time 
equation. Also for completeness, we recall that a and b are real. 

It is always a little surprising that the two exit times coincide. It also seems a little strange that 
the reflection time delay (t R ), compared to instantaneous reflection, depends upon the barrier width 
I. Is the assumption of single reflection and transmission peaks correct? This depends critically upon 
the size/width of the incoming wave packet. Before discussing further what seems a result contrary to 
the infinite multiple waves described in the previous Section, we must make the following important 
observation 

oo oo 

T = T n and R=^R n . (23) 

n— 1 n— 1 

Indeed while we have begun our analysis from the two-step calculation we could have arrived at exactly 
the same result by expanding the denominators of R and T in an infinite series. The treatment of 
R and T as single wave packets represents the limit situation in which all the infinite reflected wave 
packets overlap and similarly for the transmitted wave packets. Plane waves can in this sense be 
considered as infinitely extended wave packets and they thus satisfy automatically this coherence 
condition. The two approaches are perfectly equivalent. A single peak may break up under suitable 
conditions into multiple peaks, or equivalently, multiple peaks may coalesce under suitable conditions 
into a single peak. 



V. RESONANCE PHENOMENA 

One of the characteristic of the single wave packet situation is manifest in the expressions for R and 
T. The reflection coefficient R vanishes when sin[gZ] = 0. It follows that for values of I such that 

I = nir / q (n = 0, 1, ...) , 

we obtain complete transmission. In the figure 1, we show, for different values of V , E and m, the 
typical resonance curves for |T| (dotted lines). In principle it extends to infinite values of I. Any 
normalizing convolution integral will modify this. First, an integration over p, unless very "tight" 
about any p value will imply some averaging of this curve. For a spread in momentum Ap such 
that Aq > 7r q the resonance effect will be completely averaged out. Secondly, for an incoming wave 
packet with finite spatial spread, say Az ~ l/Ap, we can ask when the multiple peaks, described in 
the previous Section, are well separated. This occurs when the distance between two peaks is much 
greater than Az. Thus, for complete decohcrence 

1 

v g . ul At > Az ~ — , 
Ap 

where v g>III is the group velocity in region III and At the time interval between successive peaks. 
Now At = 2l/vg t u, hence multiple peaks will be clearly separated when 

21— — > Az or — - > — — = — . 

v g ,u Az v g . IU p 

For a plane wave Az = oo. Hence, for plane waves we have maximum coherence always. For any finite 
(Az), we see that coherence is lost as I increases. Decoherence implies (as proven in the previous 
Section) that the total transmission probability becomes independent of / and is given by (|17|l , 



E 

n=l ~ 1 

This value happens to coincide with \T\ (but not \T\ 2 ) at the mid resonance values where cos[gZ] = 0. 
All of this is exhibited in our numerical calculations, in the figure 1, where the exact value of the 
transmission probability is plotted against I. The tendency to a constant value as decoherence sets in 
is apparent. 
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The condition for decoherence is obviously achieved as I — > oo. It is also obtained if Az — > 0. 
However, Az — ► implies Aq — ► oo and we must be careful not to drop below the A (above barrier) 
zone. There is also a third limit in which it occurs, when Vg ' 11 — > 0. Let us now consider this last 

V 9,HI 

limit in more detail. It can be achieved in two different ways. The first is by sending q — > with 
V Q fixed, whence p — > \JV a (V a + m). The second is by keeping q fixed and sending p a (E ) and V 
simultaneously to infinity so that E a — V remains constant. The first choice is again difficult to realize 
because of the Heisenberg uncertainty principle. To see this concretely, consider a symmetric convo- 
lution function about q (Po)- As q — > so to must Aq — > since we must stay above the tunnelling 
zone. This means that automatically we must have Az — > oo. So contrary to our expectations we end 
up in the coherent state. The second choice however does indeed lead to complete decoherence since 
it can be achieved while keeping the widths of the wave packets fixed. 



VI. CONCLUSIONS 

In this paper we have considered diffusion of an incoming wave (or wave packet) with E > V a + m 
by a one-dimensional potential barrier of height V a and width I. In front and beyond the barrier, the 
potential is assumed to be zero (free space). In this study, we have employed the SPM. There is an 
inherent ambiguity in this method. Given a sum of terms, it may be applied individually to each or 
to the sum. In the former case a scries of wave packet peaks are determined while in the latter case 
only one peak is predicted. It is easy to use the SPM in both the limit of complete coherence and 
complete decoherence. It is difficult to see how to use it for intermediate situations. In these cases we 
can fall back upon pure numerical calculations or possible use cluster decompositions in which the 
wave packets are summed in finite numbers before applying the SPM. However, this possibility and 
its viability has still to be explored. 

The overall reflection and transmission amplitudes (R and T) are characterized by resonance 
oscillations and by the feature that the reflection delay time is equal to the transmission time; at 
least when (or to extent) that we can consider each a single wave packet. What we have shown in 
this paper is that, even with the Dirac equation, the barrier results can be obtained by treating the 
barrier as a two step process. This procedure involves multiple reflections at each "step" and when 
I — ► oo predicts the existence of multiple (infinite) outgoing peaks. However, by simply summing the 
individual amplitudes, one obtains exactly the standard barrier results. This has lead us to postulate 
and then confirm numerically that with increasing / the resonance curves will lose coherence and tend 
to predicted constant values. 

From an alternative, but equivalent viewpoint, it has also been noted that, whereas one cannot 
perform the limit I — > oo in the R and T amplitudes, one can expand the denominators, in a natural 
way, into an infinite scries which reproduces exactly the two step results. For any finite but sufficiently 
large I we predict the appearance of multiple peaks of which the first reflected term is simply the 
single step result characterized by "instantaneous" reflection. As I increases the exit times of the 
other peaks grow. This suggests that for finite times (or simply ignoring secondary peaks) the single 
step is equivalent to a barrier with a sufficiently large width. So not only can we claim to have shown 
that a barrier is equivalent to two steps but, at least for the first reflected wave, that a wide barrier 
is in its turn an approximation of a single step. 

We shall call upon these results in a subsequent work or works in which we consider the tunnelling 
energy zones and the Klein zone. We consider these energy zones separately because they are charac- 
terized by different physical phenomena (resonance, tunnelling, pair production). For tunnelling, we 
shall be particularly interested in the extension of the Hartman effect [9-11] from the Schrodingcr 
equation to the Dirac equation [7,8]. This would justify the renaming of the effect to the Hartman 
paradox since it implies the existence of supcr-luminal velocities. Within the Klein zone we will again 
use the two step method described in this paper. 
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APPENDIX. 



In order to obtain the expression of the transmitted probability for the step potential case, we consider 
the following incident wave packet, 



$ Itl {z,t) = / dp g(p) 



/ 1 \ 



p 

, E + m , 

V o / 



exp [i{pz — Et)] , 



(24) 



where p min = y/V (V + 2m) and g{p) is a real function with a pronounced peak about the value 
p = p chosen by construction such that 



dz |<M*,t)| = / dp g\p) 



The transmitted wave packet can then be written as 



V 



E + m 



= 1 



$„{z,t) = / dp g(p)T+(Q) 



( 1 \ 



9 



Pmin 



E-V + m 




exp[i (q z — Et)] 



T + (0) 



<1 



E-V + m 




cxp[i q z\ x 



P = Po 



dp g(p) exp 



q E 



Consequently, 



(25) 



(26) 



dz \$n(z,t)\ c \T+(0)\ 
T + (0)\ 2 



E — V q E 



o E -V + m p (E - V Q ) J p 

2{E - V ) g E E a + m 
o E a - V + m p a (E - V ) 2E Q 



dp g 2 (p) 



\T + (0)\ - 
a 



(27) 
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Fig. 1. Barrier width dependence of the transmission probability for different values of aV , aE and ma 
(h = c — 1). The transmission probability |T(p )| 2 (see dotted lines) shows the typical resonance curves for 
plane waves. In our numerical calculations, performed for asymptotic times, we have used a gaussian incoming 
wave packet of width a. The transmission probability a J dpexp[— a 2 (p — p ) 2 ] \T(p)\ 2 (see solid lines) 

exhibits the tendency to the constant value 2a a b /(a 2 + b 2 ,) as predicted by the multiple scattering analysis. 



